Draft version February 5, 2008 

Preprint typeset using 1^1^^ style cmulatcapj v. ^jllj^A, 



O 

o 



> 

vn 

(N 
O 

\o : 
o ■ 

Oh- 
I 

o 



A METHOD TO IDENTIFY THE COMPANION STARS OF TYPE lA SUPERNOVAE IN YOUNG 

SUPERNOVA REMNANTS 

Jin Ozaki^' ^ and Toshikazu Shigeyama^ 

Draft version February 5, 2008 

ABSTRACT 

We propose a method to identify the companion stars of type la supernovae (SNe la) in young 
supernova remnants (SNRs) by recognizing distinct features of absorption lines due to Fe I appearing 
in the spectrum. If a sufficient amount of Fe I remains in the ejecta, Fe I atoms moving toward us 
absorb photons by transitions from the ground state to imprint broad absorption hues exclusively with 
the blue-shifted components in the spectrum of the companion star. To investigate the time evolution 
of column depth of Fe I in the ejecta, we have performed hydrodynamical calculations for SNRs 
expanding into the uniform ambient media, taking into account collisional ionizations, excitations, 
and photo-ionizations of heavy elements. As a result, it is found that the companion star in Tycho's 
SNR will exhibit observable features in absorption lines due to Fe I at AA = 3 71.9935 nm and 385.9911 
nm if a carbon deflagration SN model l)Nomoto. Thielemann, fc Yokoilll984D is taken. However, these 
features may disappear by taking another model that emits a few times more intense ionizing photons 
from the shocked outer layers. To further explore the ionization states in the freely expanding ejecta, 
we need a reliable model to describe the structure of the outer layers. 

Subject headings: general — supernova remnants: individual(Tycho's SNR) 



1. INTRODUCTION 

Type la supernovae (SNe la), characterized by no Hq 
but strong Si lines in the spectra at the maximum bright- 
ness, are brighter than most of SNe classified into the 
other types and exhibit uniform light curves. Thus they 
are used as a standard candle to measure distances to 
remote galaxies. 

A plausible explosion model for SNe la is the accret- 
ing white dwarf (WD) model, in which a white dwarf in 
a binary system accretes material from the companion 
star, increases its mass, usually up to the Chandrasekhar 
mass limit (Mch ~ 1-4 M0), and then explodes (e.g., 
iNo moto 1982). There have b een significant progresses in 
the ac creting WD model since lHachisu. Kato. fc Nomotd 
l)1996j) introduced the stellar wind from the WD while it 
accretes materials from the companion. Their model suc- 
ceeded in sustaining a stable mass transfer in the progen- 
itor systems of SNe la. According to their model, there 
are two main evolutionary paths leading to SNe la, the 
super soft X-r ay source (SSS) chann el and th e symbiotic 
channel (Hachisu. Kato. fc Nomo tciri 996: Ha chisu et al.l 
Il999t IHachisu. Kato. fc Nomotd Il999fl . Accordingly 
their model predicts which companion stars lead to SNe 
la. 

The above evolutionary scenario for SNe la has not 
been confirmed by observations, which will require the 
identification of the companion star that should remain 
in the vicinity of the explosion site. Ruiz-Laoucntc ct al. 
l)2004tl argued that their group identified a G2IV star 
as the companion of Tycho Brahe's supernova remnant 
(SNR) by measuring the velocities and distances of stars 
in the vicinity of the center of the SNR. They concluded 
that this G2IV star was moving much faster than the 
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other neighbor stars and that the distance to this star 
seemed consistent with the distance to Tycho's SNR. Al- 
though Fuhrmann (200^ has argued that the observed 
velocity of Tycho G might correspond to the velocity of 
stars belonging to the thick disk population, the expected 
stellar mass in thick disk stars within the cone with 2.87 
arcsec radius (which corresponds to the angular distance 
from the center of Tycho's SNR to the Tycho G star), 
at 3 kpc from Earth, is only 2 x10~'^Mq, which makes 
the thick disk star alternative very unlikely. For this es- 
timate we use the density in the vicinity of Tycho's SNR 
ijDehnen fc Binnevlll998|i . 

Although the coincidence of the kinematic character- 
istics of Tycho G with its being at the position and dis- 
tance of the SNR appears significant, confirmation by 
other means would nonetheless be very useful. 

In this paper we propose a direct method to prove that 
the companion star is located inside the SN ejecta. A 
hint was dropped by observations for a star called S-M 
star discovered by Schwcizer fc M iddle ditch (1980) near 
a type la SNR 1006. IWu et all lfl993f) proved that the 
S-M star was not the companion of this supernova by in- 
vestigating features of Fe II absorption lines in the ultra- 
violet (UV) spectrum. Very broad wings were observed 
in both blue and red sides of the absorption lines. The 
line width was a few thousand km/s much larger than the 
thermal velocity of stellar atmosphere, which is thought 
to be '^10 km/s. The broad wings are likely to be formed 
by Fe II in the ejecta of SNR 1006; photons in the blue 
wing are absorbed by the matter ejected toward us, and 
those in the red wing away from us. Thus it was proved 
that the S-M star is located behind the SNR 1006 ejecta. 

When a star is inside the ejecta of SNe la, photons 
emitted from the star are absorbed only by the ejecta 
moving toward us. Hence the broad wing must be present 
only in the blue side. The absorption line with only blue 
wing enables us to identify companion stars of SNe la. 
IWu et alJ l|1993f) used the UV range to observe the S-M 
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star with the Faint Object Spectrograph on the Hub- 
ble Space Telescope. However, in addition to difhculties 
in UV observations from the ground, companion stars on 
the evolu tionary paths suggested by the above mentioned 
scenario (Hachisu. Kato. & Nomoto 1996: Ha chisu etaP 
IT999: Hachisu. Kato. fc Nomoto 1999) may not be bright 
in the UV range. Then we will focus on absorption lines 
in the visible range. Furthermore the corresponding tran- 
sitions need to be from the ground state because most Fe 
ions in the ejecta are expected to be in the ground state. 
Thus only Fe I can produce such absorption lines in SNe 
la ejecta. 

In this paper, we estimate the amount of Fe I in the 
ejecta of Tycho's SNR by taking account of collisional 
processes in non-equilibrium and ionizations by photons 
emitted from the shocked ejecta, and calculate spectra of 
a star located at the center of a SNR and discuss whether 
we can identify the feature of Fe I absorption lines in the 
spectrum of the companion star. 

2. FORMULATIONS 
2.1. Hydrodynamical Model 

We calculate the hydrodynamical evolution of a spheri- 
cally symmetric SN la embedded in the ambient medium 
with a uniform density. The initial structure of 
freely expanding SN la ejecta is t aken from W7 in 
INomoto. Thielemann, fc Yokoil l(T98^ . The mass of the 
WD at the time of the SN explosion is 1.37 Mq. Our 
model simplifies the original distribution of elements in 
W7. That is, the ejecta consist of three layers with differ- 
ent compositions; the innermost layer is consists of pure 
Fe with the mass of 0.95 Mq, the middle layer consists of 
0.15 Mq Si, and the outermost 0.27 Mq C/0 . The mass 
fraction of oxygen in this layer is 50 %. The initial kinetic 
energy is 1.3 x 10^^ ergs. The densit y rta of the ambient 
medium is assume d to be 0.41 cm~^ llHamilt on fc FesenI 
[1988) or 1.0 cm-3 ijltoh. Masai, fc No moto 1988). These 
two values refer to two different models of Tycho's SNR. 

To obtain the time evolution of physical states in the 
ejecta, we solve hydrodynamical equations with the La- 
grangian coordinate 
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where r is the specific volume, u the velocity, E the total 
energy density per unit mass, m the Lagrangian mass co- 
ordinate. The internal energy e, the mass density p, and 
the pressure p are derived from the conserved quantities 
via 
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where 7 is the ratio of specific heats and assumed to be 
5/3 for simplicity. Here we have ignored the radiative 
cooling, because the associated time-scale is significantly 
longer than the exp ansion time-scale of a young SNR like 
Tycho's SNR fe.g.. ISorokina et al J 120041. We numeri- 
cally integrate the above e quations ^ with the PPM 
(IColella fc WoodwardllT98l . 



2.2. Two Temperature Approximation 
for Shocked Plasma 

In the shocked ejecta, electrons and ions are not imme- 
diately relaxed to an equilibrium state because of their 
different masses and low densities. Two temperature ap- 
proximation is used to describe this relaxation process. 
The time scale for the electron temperature To t o catch 
up wi th the ion temperature T\ is estimated by 'SpizgJ 
l)1978(l . The time scale teq, which is of the order of 10* 
years, is longer than the age of Tycho's SNR. Then we 
must calculate the electron temperature by solving the 
equation 

m ^ T,-T, 

dt ioq(i,e)' 

2.3. Collisional Processes in the Shocked Region 

In the shocked region, the electron temperature be- 
comes higher than 10*" K. Thus collisional excitations 
and ionizations are activated and dominated. To obtain 
the population f{X^) of an element X in the ionization 
stage J, we solve rate equations written as 



dt 



= p-'fix^-')-pf{x^), 



(3) 



where P denotes the rate of ionizations out of the jth 
ion stage. The collisional ionization cross sectio n is cal - 
culated using the formula introduced in iLotd l)1967j) . 
The excitation-autoionization (E-A) is taken into ac- 
count for ions such as Fe^+, Fe^+, Fe^°+, Fe"+, Fe'^'^+, 
Fei3+. Fei*+. Fe'^+ . Ff^+ Si^+, Si"+, 0^+, and 
ijArnaud fc Rothenfluelll985(l . On the other hand, we 
can ignore recombination rates in the above equations, 
because the recombination time scale becomes of the or- 
der of 10^ years much longer than the age of Tycho's 
SNR. 

Instead of taking into account photo-ionizations in 
front of the shock, atoms immediately after the shock 
front in the ejecta are assumed to be automatically ion- 
ized to Fei°+, Si'^+, 0^+, 0"^+ to account for the detec- 
tion of Si XIII and Fe XV II lines by the XMM-Newton 
l)Decourchelle_et aT 200 1|). This assum ption has been 
made bv lltoh. Masai, fc Nomotd 1)1988(1 . The photo- 
ionization process is taken into account as a post-process 
to reduce the computational time. 

2.4. Photon Emission 

A lot of high energy photons are e mitted from ions of 
Si, Fe, O, and C i n the shocked ejecta l|Hamilton fc Fesei] 
IT988tlVink eTal]l2003|) . These photons ionize the matter 
in the pre-shocked ejecta. The number of ionizing pho- 
tons emitted from ions i per unit volume per unit solid 
angle per unit time is 
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where Ui is the number density of ions i, and C, 
the collisional excitation coefficient for the correspond- 
ing transition llKa.tnin 97(11) . 

We include ionizing photons emitted from transi- 
tions of ions listed in Hamilton fc Fcscn (1988) and 
IHamilton fc SarazinI 1(1984^) . Strong emission lines from 
carbon ion s {qf > 10"^) a re taken from iEngstrom et alJ 
Ipp) and lTunklevI llT997tl . 
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TABLE 1 

Fe column densities for SNRS at the age of 430 YEARS. 



Density of 


Fe I 


Fe II 


Fe III 


Fe IV 


Fe V 


ambient medium 


(10^5 cm" 


-2) (10^5 cm" 


-2) (10^5 cm- 


2) (IQl^cm- 


2) (10i5cm-2) 


1.0 cm"-^ 


1.2 


6.0 


9.5 


22 


18 


0.41 cm~^ 


8.9 


28 


16 


8.4 


1.9 



2.5. Radiative Transfer 

With the source term for ionizing photons discussed in 
the preceding subsection, we solve the radiative transfer 
equation for ionizing photons to estimate the ionization 
states of Fe, Si, O, and C in the pre-shocked ejecta. The 
intensity / of ionizing photons is obtained by solving the 
transfer equation, 

dl l~^i^dl ^ 

+ ^7 =^Ph-Z^^Pi,cl"cH/, (5) 

where r denotes the distance from the center, 9 the angle 
between the direction of incident photons and the radial 
direction, /i = cos 9, (Tpi_ ei is the photo-ionization cross 
section of Fe, Si, O, or C ((Verncr ct al. 1996), and 'el, i' 
represents (* — 1) times ionized element el. This equation 
is integrated along the characteristics to obtain I{r, 9) 
in th e freely expand ing ejecta as described in textbooks 
fe.g.. iMihalasI Il978^ . This equation is solved after the 
optical depth becomes smaller than ~ 10^. 
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Fig. 1. — Radial profile of X-ray surface brightnesses of Si emis- 
sion (1670-2000 eV; solid fine) and Fe XVII emission line (A = 15 A; 
dashed line) from a model with = 1.0 cm~'^ at a distance of 2.43 
kpc from the earth. 



2.6. Ion Fraction 

The number /ph of photons incident from all directions 
per unit time is obtained by integrating the intensity / 
with respect to angle. If the corresponding value for pho- 
tons that can ionize ions is denoted by /ph, x^j then 
the ion fraction f{X^) of element X in the jth ioniza- 
tion stage evolves with time according to the following 
equations, 

= f{X^-')ia^,, x.-^ /ph, X.- + n,{v<Jcoiy-') 
-f{X^){(Jpi^xifph,xi +nc{vaco\y) (6) 

where cTpj is the cross section for photo-ionizing jth 
ions of an element X and nc{vaco\Y is the collisional 
ionization rate. 

3. RESULTS 
3.1. X-ray Surface Brightness 

Comparing the calculated X-ray surface brightness 
profile with that obse rved for Tycho's SNR (e.g., 
iDecourchelle et aLll2nni|) . one can constrain the density 
in the ambient medium. Figure shows the calculated 
X-ray surface brightness of Si XIII triplet blend (in the 
range of 1670 - 2000 eV) and Fe XVII hue (at A = 15 A) 
emitted from a SNR model with — 1.0 cm~'^ located 
at a distance of 2.43 kpc from the earth. The plotted 
profile has been smoothed by a Gaussian filter of radius 
5" and binned to a 12.4" interval to mimic the observa- 
tion. The observed radius of the peak is reproduced quite 
well by this model, while the contrast between the peak 



brightness and the central brightness is slightly more pro- 
nounced in the model. This might be due to a devi- 
ation from spherical symmetry in Tycho's SNR, which 
will lower the peak brightness when plotted assuming 
spherical symmetry. 

3.2. Column densities of Fe ions 

The column densities of Fe ions thus obtained for SNRs 
at the age of Tycho's SNR are listed in Table [T] 

Comparing with those of Hamilton & Fcscn (1988), 
our results show that ions tend to populate in lower ion- 
ization stages in all calculations. There seem to be two 
differences in these two calculations. The main difference 
is in the treatment of carbon ions as a source of ionizing 
photons. Hamilton & Fcscn (1988) replaces carbon with 
oxygen in calculating the number of ionizing photons, 
while we do calculate the number of ionizing photons 
from carbon ions. The other is a difference in hydrody- 
namical models. We use the W7 model, but they use 
the C DTG7 model constructed by Wooslcv & Weave^ 
(1986). Because the CDTG7 model has a less centrally 
condensed density profile, photo-ionizing emission from 
the outer layers are more intense than in the W7 model. 
Both of th ese two difference s enhance photo-ionizations 
in Hamilt on fc Fesei] (|1988|) as compared with our cal- 
culation. 

4. SPECTRA OF THE COMPANIONS 
4.1. Radiative Transfer 
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Some photons emitted in the visible band from the 
companion star are absorbed by Fe I in the expanding 
ejecta. Suppose the companion star emits photons with 
the frequency Vs toward an observer. These photons are 
absorbed by Fe I moving with the velocity v, when 



(7) 



is satisfied as long as v is much smaller than the speed of 
light c. Here f o denotes the frequency of Fe I absorption 
line in the co-moving frame. Thus the absorbed photons 
must have a frequency higher than vq. 

The optical depth t{vs) of the pre-shocked ejecta for 
photons with the frequency is given by 



(^ahsnpcidr. 



(8) 



Here the radius of the inner shock front is denoted by R. 
This integration can be transformed into the integration 
with respect to frequency as 

r=(i-«/=*) /dr\ /dv\ , 

<7s^bsnFcldr = / O-absnpcI ( ^ ) ( ^ ) 

(9) 

where v denotes the photon frequency in the co-moving 
frame, i.e., ly — Vsi^ — v/c), crabs the absorption cross 
section, and npci the number density of Fe I. 

Since the velocity of the ejecta (about several thousand 
km/s) is much larger than the velocity of the thermal 
motion of particles in the pre-shocked ejecta (less than 
a ten km/s), applying the Sobolev approximation to the 
integration above yields, 



ne 
rricC 





gfnFci 



' dr ^ 
V dv J 



I I 
I du I 



if vo <Vs < Vq{1 + R/ct), 
otherwise. 

Here the density npci and the derivative dv / dv are eval- 
uated at the point where vq = 1^3(1 — r/ct) is satisfied. 

We investigate the transfer of photons in the U band 
absorbed by four transitions of Fe I listed in Table |2] 

Though we have assumed that the distribution of Fe is 
a function of the distance from the center, Fe may be dis- 
tributed in blobs. We will briefly discuss how these blobs 
affect the absorption lines discussed in this section. Sup- 
pose Fe is distributed in N identical blobs with spherical 
shapes. These blobs fill a fraction / of the entire volume 
of the ejecta. Then the column depth of Fe is reduced 
when 

3 



^^2/3^1/3 < 



(11) 



is satisfied. Otherwise the column depth will not be 
changed but the shape of the absorption lines will be 
changed depending on where blobs on the line of sight 
are loc ated in the ejecta. A 3-dimensional modeling of 
SNe la l|Reinecke et al.ll200^ demonstrated that a larger 
number {N ^ 10) of bubbles is preferred to accommodate 
the observed amount of ^^Ni (Fe) and explosion energy. 
Thus a large N expected from this 3-dimensional mod- 
eling is unlikely to reduce the column depth. 

dumpiness in the ejecta was suggested by recent ob- 
servations for SN 2004dt fWane et al. 200|). Spectropo- 
larimetry of this supernova more than a week before its 
optical maximum reveals that the Si layer is substan- 
tially deformed from spherical symmetry while the same 
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Fig. 2. — (a) Spectra of a G2IV star located at the center of 
a SNR at the age of Tycho's SNR. The long dashed line shows 
the spectrum for = 1.0 cm~'^ and the short dashed line for 
ria = 0.41 cm~^. The solid line exhibits the original spectrum of 
the star, (b) the same as (a) but for a star embedded in clumpy 
ejecta described in section 4.1. 



observation shows little polarization signature in O lines. 
Since Si is not a major source of ionizing photons, the in- 
fluence to the Fe I lines is insignificant. From the theoret- 
ical point of view, howev er, the clumpy distrib ution of Fe 
I can not be ignored fe.g. lReinecke et an i2002'l. Then we 
will calculate the spectra from a clumpy distribution of 
Fe I in the next section though the distribution is rather 
artificial. To mimic clumpy structures especially in the 
inner region of the ejecta expected from results of 3-D 
calculations, the density of the innermost ejecta with ve- 
locity <2,000 km/s is reduced by a factor of two and the 
reduced mass is added to the adjacent outer zones with 
the velocity range of 2,000-3,000 km/s by increasing the 
densities thereof by a constant factor. Here we assume 
that the ion fractions are not affected by the dumpiness. 

4.2. Spectra 

The spectra of a G2IV star located at the center of 
SNRs embedded in the ambient media with two different 
densities are calculated and shown in Figure |2 (a). This 
star has t he same spectra l type as the candidate star 
Ruiz-Lapu ente et alJ l)2004D found. 
Figure |2Ib) shows the spectra of the same star in the 



ejecta in which the distribution of Fe I is affected by 
the clumping discussed in the previous subsection. It 
is clear that the clumpy distribution of Fe I changes the 
shapes of absorption lines: the lower velocity components 
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TABLE 2 
Equivalent widths of Fe I lines 



Fe I Line 




Equivalent widths 


Wave Length 


log(9/) 


(nm) 


(nm) 




ria =1.0 cm 


ria = 0.41 cm~^ 


367.9913 


-1.599 


0.363 


1.95 


371.9935 


- 0.431 


2.51 


4.88 


382.4444 


-1.362 


0.633 


2.86 


385.9911 


-0.710 


1.94 


4.59 



of the hnes are weaken because of the reduced density in 
the corresponding layers while the addition of mass to 
the outer layer hardly affects the line shape because a 
substantial amount of Fe is already ionized there. The 
equivalent width of each line exhibited in Figure |2] (a) 
is listed in Table |21 This particular dumpiness in Fe I 
distribution reduces the equivalent width by about 30 %. 
The spectra exhibit strong blue wings both for tenuous 
and dense ambient media. Results of our calculations 
indicate that Tycho's SNR has a sufhcicnt amount of Fc 
I to identify the companion, even though a companion 
star as late as G2 has a very complicated spectrum with 
many absorption features. 

Even when we find a star that exhibits the absorp- 
tion feature discussed in this paper, there is a chance 
that a star other than the companion star happens to 
be inside the ejecta. A star existing in the vicinity of 
a SN la gets a fraction of the explosion energy and is 
accelerated. If the companion star similar to Tycho G 
star with the mass of Mc = 1.2 Mq and the radius of 
Rc — 2 Rq is located at the distance of i?o = 5.5 Rq 
from the progenitor, the size of the Roche lobe is com- 
parable to the size of the companion. Thus the star in 
this situation will get the maximum velocity after the 
explosion. Then the orbital velocity before the explosion 
is ~ 160 km/s. In a 3-D hydrodynam ical simulation of 
SNe la bv [Marietta. Burrows, fc Frvx ell ( 2000), they ob- 
tain a plausible kick velocity i; 50 km/s. Therefore, 
including the orbital velocity, the velocity of the com- 
panion star becomes Vc ~ 170 km/s. If the companion 
star has been moving away from the explosion site at 
that speed, the companion star is now at the distance of 
<0.08 pc from the center. The stellar mass density in the 
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neighborhood of Tycho's SNR estimat ed from a model 
of the Galaxy ijDehnen fc Binnevlll998|) is less than 0.03 
AIq/pc~^. Thus the expected stellar mass inside this 
volume is only ^ 1O~^M0. Therefore if we find a star 
showing absorption lines with broad blue wings in the 
spectrum, it is likely that the star is the companion star. 

5. CONCLUSIONS AND DISCUSSIONS 

As a consequence, we have demonstrated that there 
exhibit very deep absorption lines with unique shapes 
in the spectrum of the companion star located at the 
center of a young SNR such as Tycho's SNR. There are, 
however, a few factors that might reduce or even erase 
these distinct absorption features. 

First, the number of Fe I in the freely expanding ejecta 
is very sensitive to the number of ionizing photons emit- 
ted from the shocked ejecta. An increase in the number 
of ionizing photons by a few factors might decrease the 
number of Fe I by a few orders of magnitude or more. 
Since the main source of ionizing photons is O in the 
outer ejecta, the distributions of O and density in the 
outer ejecta need to be known precisely. Unfortunately, 
these regions in W7 have s ome problem s to reproduce 
the observed optical spectra |Bra nch| 1998l and references 
therein). Due to this uncertainty in the explosion model, 
it is not conclusive if the companion star of Tycho's SN 
will exhibit unique Fe I absorption features discussed in 
this paper. Nevertheless, it is true that every SN la has a 
period during which the companion star has the distinct 
absorption features discussed above because the ejecta 
become cool enough to have plenty of Fe I for a time 
after the optical brightening. 

Second, it is assumed in our calculations that ions in 
the shocked region are ionized to Fe^°+, Si^+, 0^+, C'*'*' 
immediately after the shock passes fo l lowing the proce- 
dure taken bv lltoh. Masai, fc Nomotcl l)1988(l . Since ions 
in lower ionization stages are a strong source of ionizing 
photons, ionization may be more advanced in real young 
SNRs. 
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